Formation and local symmetry of Holstein polaron in t-J model 



(N 

O 
(N 

o 

Q: 



Han Ma\ T. K. Lee^, Yan Cheni^^ 
^Department of Physics, State Key Laboratory of Surface Physics and 
Laboratory of Advanced Materials, Fudan University, Shanghai 200433, China 

'^Institute of Physics, Academia Sinica, NanKang, Taipei 11529, Taiwan 
'^Department of Physics and Center of Theoretical and Computational Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong, China 
(Dated: December 5, 2012) 

The formation and local symmetry of spin-lattice polaron has been investigated semiclassically in 
the planar Holstein t-J-like models within the exact diagonalization method. Due to the interplay 
of strong correlations and electron-lattice interaction, the doped hole may either move freely or 
lead to the localized spin-lattice distortion and form a Holstein polaron. The formation of polaron 
breaks the translational symmetry by suppression of antiferromagnetic correlations and inducement 
of ferromagnetic correlations locally. Moreover, the breaking of local rotational symmetry around 
the polaron has been shown. The ground state is generically a parity singlet and the first excited 
state maybe a parity doublet. Further consequences of the density of states spectra for comparison 
with future STM experiments are discussed. 

PACS numbers; 74.78.Fk, 71.10.Fd, 74.20.Fg 



k 



I 

O 

o 



>: 
o: 

^ ' 



X: 



Doping a Mott insulator has been regarded as the main 
physics in high Tc cuprate superconductors 1]. Single 
hole in the two-dimensional t-J model may form a fer- 
romagnetic (FM) spin polaron for sufficiently small ex- 
change coupling J 0. For finite J, the distortion of 
the antiferromagnetic (AFM) background decays away 
from the hole and the competition between magnetic 
correlation and kinetic energy may result in the ground 
state with a spin polaron structure. This problem has 
been extensively investigated numerically [3|-|7[ • To make 
a better comparison with angle-resolved photoemission 
spectroscopy (ARPES) experimental data, one need to 
include the nearest-neighboring (NN) and next-nearest- 
neighboring (NNN) hopping t' and t" terms (sl-fiot. On 
the other hand, the role of electron-phonon (el-ph) cou- 
pling has gained much interest recently. One reason is 
that the ARPES data in doped metallic cuprates which 
showed the broadening of spectral lines at a certain mo- 
mentum revealed the band dispersion renormalized by 
e^-ph interaction [ll|, [l^l- In addition, the interaction 
also shifts the energy of the states. 

The doped charge carriers in the presence of both 
strong electronic correlations and electron-phonon inter- 
actions may lead to the formation of spin-lattice polaron. 
In particular, the AFM exchange interaction allows for 
spin flips leading to coherent hole motion at the bottom 
of the band and forms a spin polaron. In the presence 
of strong eZ-ph coupling, both the spin and lattice de- 
grees of freedom become entangled and the spin polaron 
may transform into a spin-lattice polaron. The forma- 
tion of this composite polaron may affect both the spin 
and lattice degrees of freedom locally. Recent ARPES 
experiments in undoped cuprates were interpreted in 
terms of strong eZ-ph coupling giving rise to localized po- 
laron [111, [13 ■ The possibility of self-localization of holes 
in lightly doped cuprates has been studied [l^. It has 
been found that the effect of eZ-ph interaction on spin 



polaron is strongly enhanced as compared to polaron in 
uncorrelated systems p^ - [l9| . 



Recent atomically resolved scanning tunneling mi- 
croscopy (STM) studies [l^, [2l| on strongly underdoped 
cuprates revealed a surprisingly pattern with the square 
symmetry of the lattice broken on a local scale. The 
origin of this broken local symmetry was attributed to 
the dopant impurity effect '22]. The broken symmetry 
states were shown to appear in the case of a hole confined 
to a cluster of sites centered at a impurity. Meanwhile 
the introduction of eZ-ph interaction to the t-J model 
may stabilize the half-doped stripes [1^ . In the presence 
of strong eZ-ph coupling, the variations of hopping inte- 
gral and spin-spin correlation around the impurity may 
become more remarkable and the composite spin-lattice 
polaron may show up. 



In this paper, we shall discuss the formation of spin- 
lattice polaron and its relevance of local symmetry by 
investigating the Holstein t- J-like models with the exact 
diagonalization method. Due to the interplay of com- 
peting electronic correlations and eZ-ph interactions, the 
doped hole may either move through lattice freely or fa- 
vor the composite spin-lattice polaron. The formation 
of Holstein polaron breaks the translational symmetry 
by suppression of AFM correlations and inducement of 
FM correlations locally. Moreover, the breaking of lo- 
cal rotational symmetry around the polaron center has 
been shown. The ground state is generically a parity sin- 
glet and the first excited state maybe a parity doublet. 
Further consequences of the tunneling spectra for com- 
parison with future STM experiments will be discussed. 



The two-dimensional single-band Holstein t-t'-J model 
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FIG. 1: The total energy (a) and kinetic energy (b) as a func- 
tion of el-ph coupling constant g. The solid line corresponds 
to the delocalized state while the dashed line denotes the lo- 
calized Holstein polaron state. The inset shows a measure of 
polaron size (d) as a function of g. We choose t' = —0.1 and 
J = 0.3. 



in the adiabatic limit is defined by the Hamiltonian, 
^ = E (^Ic,- + h.c.) - i' E (clcj- + h.c.) 

{i,j) i i 

where c| ^ is an electron creation operator with spin a at 
site i with a constraint of no double electron occupation 
due to strong electron correlations, Si is a spin opera- 
tor for electron at site «, and refer to NN 
and NNN sites i and j. The first three terms in Eq. 1 
represent the usual t-t'-J model Hamiltonian. The fourth 
term denotes the e/-ph interaction. In the adiabatic limit, 
displacements Ui have been treated classically and deter- 
mined by Ui — g/K{n'-). The hole density operator is 
defined a.s — l — rii — 1 — cj ^Ci^cr- The last term is 
lattice elastic energy with elastic force constant K . The 
energy spectrum and eigenstates can be obtained through 
exact diagonalization method. The symmetry of the low- 
est energy state is sensitive to boundary condition and 
parameters. In our calculations, we choose t = K ^ \^ 
suitable for the cuprates, and adopt the periodic bound- 
ary condition. 

Employed with the exact diagonalization method for 
a finite 16-site square cluster, the low-lying electronic 
states have been calculated as a function of e/-ph cou- 
pling constant g. Due to the competition between elec- 
tronic correlations and electron-lattice interactions, the 



ground state has certain limiting cases. When g is very 
small, the ground state must correspond to a delocal- 
ized state, that is to say, doped hole may move through 
lattice freely so that the average occupation number of 
holes at each site is uniform. The presence of strong 
eZ-lattice interaction may clearly favor the localized hole 
state, which is a sliding periodic polaron lattice. Since 
we focus on the local properties of polaron, we use a very 
weak impurity to break the translational symmetry and 
pin down the sliding polaron lattice so that the doped 
single hole may stay around certain site. At g 3> 1 fimit, 
the doped hole tends to be localized and results in large 
lattice distortion around that site. The corresponding 
ground state is a localized state with polaron formation. 
The evolution of total energy as well as kinetic energy 
of the system as a function of g is depicted in Fig. 1(a) 
and Fig. 1(b), respectively. The critical e/-ph coupling 
f^c ~ 1-3 can be straightforwardly determined by distin- 
guishing these two distinct states. Moreover we use the 
quantity {d) to measure the size of polaron qualitatively, 
(d) = ^ ■ ri ■ n'^{ri), where is the distance between site 
i and polaron center. It is obvious that the larger (d) 
value may correspond to a larger polaron size. As shown 
in the inset of Fig. 1, as 5 increases, doped holes tend 
to concentrate at the polaron center so that the quantity 
{d) may be significantly suppressed. It is worth to men- 
tion that the critical point gc is rather clearly revealed 
in the behaviors of {d) and kinetic energy than that of 
the total energy. Comparing with the previous study 
on the dopant impurity effect 22;], the strong variations 
of hopping integral and spin-spin correlation around the 
polaron may lead to the formation of local lattice distor- 
tion and the appearance of tightly bounded spin-lattice 
polaron due to strong el-ph coupling [l3j . 

Next we examine the effect of NNN hopping integral t'. 
As we know, the t' term plays an important role in un- 
derstanding the superconducting correlations in cuprate 
superconductors [2J-l26j. In particular, the positive t' 
case corresponds to the electron-doped system while the 
negative t' case corresponds to the hole-doped system. In 
the following, we systematically study the dependence of 
gc as a function of t' for two distinct t-J like models. Due 
to the competing tendency between the polaron forma- 
tion with strong el-ph interaction and the itinerant elec- 
trons with large kinetic energy, we naturally expect that 
the formation of localized polaron may require stronger 
e/-ph interaction to compensate larger kinetic energy \t'\ 
term. For the Holstein t-t'-J model, it is indeed the case 
as shown in Fig. 2(a) where gc goes up as \t'\ increases. 
Another feature shown in Fig. 2(a) is the asymmetrical 
behavior between t' < and t' > regions. An intu- 
itive physical understanding of such effect can be given 
as follows: for positive i', the t' and t terms in kinetic 
energy match quite well, so that the positive t' term may 
enhance the kinetic energy more effective than the neg- 
ative t' term, hence require stronger el-ph interaction gc 
to realize the localization of Holstein polaron. 

As suggested in recent studies, the t-t'-J- J' model may 
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FIG. 2: The critical el-ph coupling Qc as a function of t' for 
Holstein t-t'-J model (a) and t-t'-J-J' model (b). We choose 
J=0.3 and J'/ J = {t'/tf. 



be suitable for iron-based superconductors [13, [13 ■ In 
such systems, the appropriate range of t' is much larger 
than that of cuprates and the NNN superexchange cou- 
pling is expressed as J'/ J = {t' /t)-^. In the case of weak 
t' term, the results are similar to that of the Holstein t-t'- 
J model. The competing nature between J' term and J 
term may lead to the effect of geometrically frustration. 
In such systems, the pairwise J' interaction does not co- 
incide with the square lattice geometry, which may sup- 
press the NN AFM correlation functions. For stronger t' 
term as well as J' term, the suppression may become so 
strong that the local FM correlation may emerge around 
the polaron center. And it may result in the enhance- 
ment of kinetic energy and the reduction of the critical 
e/-ph interaction gc, which makes the formation of po- 
laron easier. These relationships are clearly exhibited in 
Fig. 2(b) where gc decreases as t' becomes quite large. 

Similar to the dopant impurity case (22j . the presence 
of a localized spin-lattice polaron may affect not only the 
local charge and local spin distributions but also the sym- 
metry of the ground state wavefunction. In the following, 
we adopt the same parity symmetry to characterize dif- 
ferent quantum states 22]. In particular, we focus on the 
reflection symmetries of a two-dimensional square lattice 
with respect to x- and y-axes passing through the center 
of polaron (P^^ and Py respectively) and on the parity 
Since [Px{y),H] = 0, we may classify states ac- 
cording to the quantum numbers of Px, Py We denote 
the state with {P^ = +1, Py = +1) as state (++), doubly 
degenerate state (-1-1, —1), and (—1,-1-1) as states (H — ) 

and ( — h), and (—1, —1) as state ( ). As we know, in 

the absence of el-ph interactions, for a 16- (4x4) site clus- 
ter with periodic boundary condition, the ground state 
of a single hole in the t-t'-J model has a four- fold symme- 
try for t' < , which can be represented by their parity 
symmetry {++),{+-),{-+),{ )■ 

As clearly shown in Fig. 3, the ground states have 4- 
fold degeneracies for weak el-ph interaction g < gc ^ 1-3, 
For g becomes slightly larger than g^ such degeneracies 
are broken. We note that the four low-lying states with 
different reflection symmetries are quite close in energy. 
The state with (H — h) symmetry always has the lowest en- 
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FIG. 3: (Color online) Evolution of four low-lying energy lev- 
els and their corresponding parity symmetries as a function 
of el-ph coupling g. We choose t' — —0.1 and J = 0.3. 



ergy, the doubly degenerated states with parity symme- 
try (H — ) and ( — h) correspond to the first excited states, 

while the state with ( ) symmetry is the second excited 

state. As g exceeds 1.5, the state (-I--I-) originally at much 
higher energy may drop down significantly and cross the 

the singly degenerated ( ) state. For g > 1.7, both two 

lowest energy states have the (-I--I-) symmetry. According 
to our numerical results, the ground state of spin-lattice 
polaron is quite robust in the (H — h) symmetry. 

To explore further the interplay between spin and 
charge degrees of freedom, we study the spatial distribu- 
tion of spin-spin correlation function and hopping inte- 
gral around the polaron. As depicted in Table I, we show 
the expectation value of spin-spin correlation function 
(SiSj) as well as hopping integral (cjcj) for six distinct 
bonds illustrated in Fig. 4(a). The interplay of strong 
correlation and el-ph interaction may lead to hole local- 
ization and result in a remarkable hopping integral and a 
weak FM spin-spin correlation function around the hole. 
In particular, the AFM correlation along bond 1 is com- 
pletely suppressed, and a weak FM correlation emerges. 
The spin-spin correlation function recovers quickly to the 
value -0.34 for farther bonds. Meanwhile the hopping in- 
tegral along bond 1 is much stronger than that of the rest 
bonds and the hopping integral vanishes quickly close to 
the lattice boundary. It is rather clear that the induce- 
ment of local FM correlations around the polaron may 
help the hole moving around more efficiently, then maxi- 
mizing the local kinetic energy. Hence the mutual coop- 
erative effect between spin and lattice degrees of freedom 
has been clearly revealed. For Holstein t-J model in in- 
finite dimensions, the interplay between the formation of 
a lattice and magnetic polaron in the case of a sin gle h ole 
in the AFM background has been studied before [20|. It 
shows that the presence of AFM correlations favors the 
formation of the lattice polaron at lower values of the el- 
ph coupling. Our numerical calculations agree well with 
their results. 

Furthermore, we investigate the parameter dependence 
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TABLE I: The spin-spin correlation functions and hopping 
integrals for various bonds in a 16-site cluster with periodic 
boundary condition at t' — —0.1, J = 0.3 and g = 2.0. The 
bond indices are labeled in Fig. 4(a). 



bond index 


1 


2 


3 


4 


5 


6 




-0.034 


0.335 


0.342 


0.337 


0.346 


0.347 




0.405 


0.018 


0.021 


0.002 


0.002 


0.000 



of the spin-spin correlation function for both t-t'-J model 
and t-t'-J-J' model. For these two models, the spin- 
spin correlation function on bond 1 and bond 2 shows 
quite different behaviors. For t-t'-J model, there is al- 
ways large AFM correlation on bond 2, while weak AFM 
or FM correlation shows up on bond 1. This strong sup- 
pression of AFM correlation locally around the center of 
polaron is due to the dramatic lattice distortion or the 
polaron formation. However, such strong suppression of 
AFM correlation on bond 2 may be significantly modi- 
fied in t-t'-J- J' model for large t' case. In such case, the 
introduction of large t' may lead to the enhancement of 
kinetic energy and the strong frustration effect due to J' 
term may greatly suppress the AFM correlation on bond 
2. Meanwhile the size of the spin-lattice polaron may be 
enlarged. 

To understand the nature of states with different parity 
symmetries, the spatial distribution of Sz has been calcu- 
lated and illustrated in Fig. 4. For delocalized states, as 
shown in Fig. 4(b)-4(e), states with four distinct parity 
symmetry (++), (+— ), ( — h), ( ) for t-t'-J model dis- 
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FIG. 4: The spatial distribution of Sz for t-t'-J model. De- 
localized states are shown in figure (b)-(e) for different parity 
symmetries. Localized states are shown in figure (f)-(i). The 
diameter of each circle is proportional to the value of local 
\Sz\- White (Black) circle represents up (down) spin. The 
star symbol labels the center of polaron. 



play quite different patterns. In Fig. 4(b), for state with 
(-I--I-) symmetry, weak up spin appears at the polaron 
center while weak AFM correlation shows up between 
the polaron center and its four nearest neighbors. In 
Fig. 4(c)-(d), for states with (-1 — ) and ( — h) symmetry, 
the asymmetry between two neighboring sites of polaron 
center along x and y directions is clearly shown. For 
localized states, the spatial distribution of Sz patterns 
are shown in Fig. 4(f)-4(i), respectively, corresponding 
to different symmetry states. For (-I--I-) symmetry state, 
weak down spin appears at the polaron center while weak 
FM correlation shows up between the polaron center and 
its four nearest neighbors. As we know, the doped hole 
may concentrate at the polaron center where the weak 
FM correlation may emerge. It is obvious that this (-I--I-) 
symmetry state is energetically favorable to the spin- 
lattice polaron formation. In Fig. 4(g)-(h), the asym- 
metry between two neighboring sites of polaron center 
along X and y directions becomes more significant than 
the delocalized states. Of all four parity symmetry states, 
this (++) symmetry state has the lowest energy. In ad- 
dition, we checked the spatial distribution of Sz in t-t'-J- 
J' model and similar results are obtained. In both cases, 
the ground state of spin-lattice polaron prefers the (++) 
parity symmetry. 

In order to check the existence of spin-lattice polaron, 
integrated differential conductance as a function of cutoff 
voltage is calculated to be compared with future STM ex- 
periments. Following the STM tunneling theory jso! ]. we 
write the integrated current at r up to a positive voltage 
V as 

/(r,c.) « Y.\^m\al,\rP^'^)\''e{uj-E^ + E^'^). (2) 

(T,m 

where a\ ^ is the electron creation operator with spin 
a at site r, |m) are eigenstates of the half-filled system 
with energy E^m w = eV , and 6* is a step function. The 
\^^^) denotes the single-hole eigenstate with energy E^^^ . 
In the following, we show the I-lo curve on various sites 
and for different eZ-ph coupling g. For convenience, we 
shift the origin of x-axis by wo = w -I- uiex where uiex = 
Eq — El^ and Eq, Eq'^ correspond to the ground state 
of half-filled system and single-hole system, respectively. 
Here we consider the contributions from low- lying energy 
states of the single-hole system. 

As we know already, the doped hole tends to stay 
around the polaron center and the hole density at the po- 
laron center is more significant than that of on other sites 
by increasing the eZ-ph interaction g. Thus we expect the 
integrated current at the polaron center would become 
larger when g goes up. This result is clearly shown in 
the left panel of Fig. 5. In the case of impurity doped 
cuprates, the conductance pattern is anisotropic as the 
tip of a tunneling microscope scans above the Cu-O-Cu 
bonds along the x{y)-axes. This anisotropy is quite pro- 
nounced at voltage around w ^ J In the present case, due 
to the formation of spin-lattice polaron, the ground state 
of one-hole system has (H — h) symmetry, the first excited 
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FIG. 5: The local integrated current up to voltage V for dif- 
ferent sites at various e-ph coupling constant <; as a function 
of oj. The left panel shows the results for signal at the center 
of polaron for different g values. The right panel depicts the 
results for signal at the two NN sites respectively along x and 
y directions. 



state corresponds to the doubly degenerated states (H — ) 
and ( — h). The presence of quadrupole interaction of two 
single-hole states or by other couplings may not change 
the symmetry of ground state but may destroy the two- 
fold degeneracy of first excited state. We consider the 
system to be in one of the two degenerate states, say in 
the state of (+— ). As depicted in the right panel of Fig. 
5, below certain cutoff-energy O.llJ), the integrated 
current shows four-fold rotational symmetry for the state 
with symmetry. For the higher cutoff-energy, the 

first excited state with (-1 — ) symmetry may play impor- 
tant role in the local symmetry breaking for the signals 
along X and y- directions. 

Our numerical results can be qualitatively understood 
in terms of the renormalized mean- field method for t-J 
model. Due to the presence of spatial inhomogeneous 
polaron, the renormalized factor gt and gs become site- 
dependent. According to previous study [3l|, the renor- 
malized factor of kinetic energy on the bond connecting 

As 



two sites i and j can be expressed as g^. 
we know, in the case of strong eZ-ph interaction, the hole 



density is highly localized at the center of Holstein po- 
laron. It is obvious that the bonds with appreciable hop- 
ping integral show up only around the polaron and may 
result in the suppression of AFM correlation locally. For 
the delocalized states, gl'' is simply a constant and so 
is the AFM correlation function. The above qualitative 
analysis agrees reasonable well with our numerical results 
shown in Fig. 3 and Table. I. 

In summary, we study both the formation and local 
symmetry of spin-lattice polaron semiclassically in the 
planar Holstein t-J-\ike models within the exact diago- 
nalization method. Due to the interplay of competing 
interactions among electronic correlations and e/-phonon 
interactions, the doped hole may either move freely or 
lead to the localized spin-lattice distortion and form a 
Holstein polaron. Since the formation of polaron breaks 
the translational symmetry, we use the parity symme- 
try with respect to the polaron center to characterize 
the localized states. The presence of spin-lattice polaron 
may suppress the AFM correlations and induce the FM 
correlations locally around the polaron. This effect may 
further stabilize the spin-lattice polaron. Moreover, this 
affect may lead to a strong localized state with par- 
ity symmetry as ground state. Moreover, the breaking of 
local rotational symmetry around the polaron has been 
shown for certain voltage cutoff. The present investiga- 
tion on the polaron formation and local symmetry may 
provide useful information for future STM experimental 
tests. 
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